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ABSTRACT 

We report a measurement of anisotropy in the cosmic microwave background radiation (CMBR) on 
7' — 22' scales. Observations of 36 fields near the North Celestial Pole (NCP) were made at 31.7 and 
14.5 GHz, using the 5.5-meter and 40-meter telescopes at the Owens Valley Radio Observatory (OVRO) 
from 1993 to 1996. Multi-epoch VLA observations at 8.5 and 15 GHz allow removal of discrete source 
contamination. After point-source subtraction, we detect significant structure, which we identify with 
emission from a combination of a steep-spectrum foreground and the CMBR. The foreground component 
is found to correlate with IRAS 100 /urn dust emission. Lack of Ha emission near the NCP suggests 
that this foreground is either high-temperature thermal bremsstrahlung (T e 10 6 K), flat-spectrum 
synchrotron or an exotic component of dust emission. On the basis of low-frequency maps of the NCP, 
we can restrict the spectral index of the foreground to f3 > —2.2. Although the foreground signal 
dominates at 14.5 GHz, the extracted CMBR component contributes 88% of the variance at 31.7 GHz, 
yielding an rms fluctuation amplitude of 82+g 2 [ 1 /iK, including 4.3% calibration uncertainty and 12% 



sample variance (68% confidence). In terms of the angular power spectrum, Ci - 
a range of multipoles / = 361 — 756, the detected broadband amplitude is 5Ti B 



+8.6 



a J" | ), averaged over 

[^+1)6^/27^/2 = 

/iK. This measurement, when combined with small angular-scale upper limits obtained at the 



59 

OVRO, indicates that the CMBR angular power spectrum decreases between I 
is consistent with flat cosmological models. 

Subject headings: cosmic microwave background — cosmology: observations 



600 and I ~ 2000 and 



1. INTRODUCTION 

In standard cosmological scenarios, baryonic matter de- 
couples from the radiation field at Zd ec — 1100. Thus, 
the horizon scale at decoupling separates scales of im- 
portance for structure formation from scales which probe 
only the primordial spectrum of perturbations. For 
z 3> 1, the Hubble radius subtends an angle A0h w 

0°.87rJo /2 (z/1100)- 1 / 2 , so that for ft = L angles £ 1° 
correspond to physical scales which were outside the hori- 
zon at decoupling (see, e.g., Kolb fc Turner 1990 ). On 



scales <; 1°, anisotropies are directly linked to causal physi- 
cal processes in the early universe. Specifically, they record 
the amplitude and phase of acoustic oscillations in the 
baryon-photon plasma and as such provide a direct mea- 
sure of the cosmological parameters which govern the early 
universe. 

In this paper, we present results from an experiment de- 
signed to detect CMBR anisotropy on 7' — 22' scales (I ~ 
600). This "RING5M" experiment is the most recent in 
a series of filled-aperture, ground-based anisotropy exper- 
iments at the Owens Valley Radio Observatory (OVRO). 
Previous experiments at the OV RO include the NCP ex- 
periment (|Readhead et al. 1989]) , in which the OVRO 40- 



meter telescope was used to place a 95% confidence upper 
limit of AT/T < 1.7 x KT 5 on power at 2 ' scales, and the 
RING40M experiment ( [Myers et al. 1993; ), which resulted 
in a detection of anisotropy, attributed to foreground con- 
tamination, of 2.3 x 10" 5 < AT/T < 4.5 x 10~ 5 (95% 
confidence) at the same resolution. 

Sin ce the COBE dete ction of a CMBR quadrupole in 
1992 (3moot et al. 1992), a large number of experiments 
have reported detections of anisotropy on > 1° scales (see 



Hancock et al. 1998, Bond 1996 for recent reviews). The 



RING5M is one of the few experiments to probe the re- 
gion of Z-space between these experiments and the high-/ 
range of the ear lier OVRO work (recent results from the 
CAT telescope ( [Scott et al. 1996 ) provide the only other 
detection on comparable scales). 

Section of this paper provides an overview of the 
RING5M experiment, while §|| and §|| review the data ac- 
quisition and relevant receiver characteristics in greater 
detail. Sections §|| and §^| describe the calibration of both 
telescopes used in the RING5M experiment; data selec- 
tion and editing are described in Results are presented 
m §0"§§l m which we also explore the data for internal 
consistency and describe tests for possible sources of sys- 
tematic error. The results of an 8.5 GHz VLA survey of 
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the RING5M fields and subsequent multi-frequency moni- 
toring of point sources for subtraction from the anisotropy 
data are presented in § |Io| , while contamination by Galactic 
foregrounds is discussed in Sections §|l^-§|l5| describe 
the detected CMBR anisotropy. The paper concludes with 
a discussion of the significance of our result for cosmolog- 
ical models. 

2. THE EXPERIMENT 

Observations were conducted at the OVRO from 1993 
to 1996, using HEMT- amplified radiometers on the 
OVRO 5.5-meter and 40-meter telescopes, fully-steerable 
parabolic reflectors with Cassegrain focus (5.5-meter) and 
primary focus (40-meter). At 31.7 GHz, the resolution 
of the 5.5-meter telescope is 7'.37 (FWHM). At 14.5 GHz, 
the OVRO 40-meter telescope provides a second frequency 
channel, for spectral discrimination of foregrounds. Since 
full illumination of the 40-meter telescope at 14.5 GHz 
would produce a ~ 2' beam, the feeds were designed to 
illuminate an ~ 11 m patch on the dish surface, provid- 
ing a good match to the the 5.5-meter beam at 31.7 GHz 
(see Table 0). The 14.5 GHz receiver is mounted in an 
off-axis configuration to minimize shadowing by the prime 
focus cage as well as scattering from the prime focus sup- 
port legs, thereby reducing the effect of ground spillover. 
This makes the beamwidth a strong function of focus po- 
sition and hence of zenith angle; the beam can vary by as 
much as an arcminute over the full zenith angle range. As 
a result, all observations at 14.5 GHz were restricted to 
lie within ±5° of the observing zenith angle for the Ring, 
■^Ring — 50° (see below). 

The 5.5-meter telescope is illuminated from the 
Cassegrain focus, so that the largest sidelobes of the feed 
illumination pattern see the sky instead of the ground, and 
the total contribution of the ground to the system temper- 
ature at 31.7 GHz is 6 K. Although the 14.5 GHz receiver 
is located at the prime focus of the 40-meter telescope, 
the under-illumination of the 40-meter dish results in a 
primary antenna pattern having its first sidelobes directed 
at the sky, and the ground contributes < 1 K to the system 
temperature. 

On both telescopes, the receiver input is continuously 
switched at 500 Hz between two feed horns separated by 
~ 22' on the sky; the recorded signal is the difference 
between successive millisecond integrations from alternate 
feeds. On the 5.5-meter telescope, the columns of air seen 
by the two feeds overlap in area by > 10% to 400 m, while 
on the 40-meter, the columns overlap by > 10% to 1 km, 
and the fast ("Dicke") switching freezes out atmospheric 
fluctuations which occur in both beams simultaneously, at 
the same time suppressing 1// noise from receiver com- 
ponents. During each measurement, azimuthal nodding of 
the telescopes between two symmetric positions, offset by 
the 22' beamthrow, provides an additional level of spatial 
switching, removing constant offsets or linear temperature 
gradients from the sky or ground (see This "double 
switching" technique has been used successfully in both 
previous OVRO anisotropy experiments. 

In the RING5M experiment, we observe 36 fields spaced 
by the 22' beamthrow in a ring around the North Celes- 
tial Pole (NCP). Field positions are given in Table § In 
order to suppress variations in the observed differential 
ground temperature introduced by telescope motion, fields 



are observed only within ±5° (±20 m ) of upper culmina- 
tion (transit, at ^Ri ng — 50°). Near transit, the separation 
of neighboring RING5M fields is approximately azimuthal, 
so that if we denote the temperature in each field by Tj, 
the quantity which results from the double switching is 
given by 

&T i =T i -±(T i - 1 +T i+1 ). (1) 

The effective beam pattern produced by the switching is 
shown in Figure ^. 
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Fig. 1. — Measured double-switched beam pattern at 31.7 GHz 
(top panel) and 14.5 GHz (bottom), obtained by scanning across 
point sources. Contours are —0.5 to —0.1 in steps of 0.1, —0.05, 
-0.03, -0.02, 0.02, 0.03, 0.05 and 0.1 to 1.0 in steps of 0.1. For 
observations of the Ring, the symmetric negative beams correspond 
to the adjacent RING5M fields. 

The interlocked geometry of the RING5M fields not only 
serves to eliminate systematic differences between obser- 
vations of different fields, but also provides a boundary 
condition ATi — which can be checked for residual 
systematics. Moreover, point sources in a RING5M field 
can be identified by the characteristic minus-plus-minus 
signature they produce in successive fields. 

3. OBSERVATIONS 

In all of the observations presented here, the telescope 
alternates the beams (referred to as the ANT and REF 
beams) on source by slewing in azimuth by an amount 
±A0 equal to the separation of the feed horns. The tele- 
scope integrates for equal times r s in each of four suc- 
cessive configurations, referred to as the A, B, C and D 
fluxes, shown schematically in Figure |^. In combination 
with the fast differencing between the feeds (ANT — REF), 
this procedure, known as a "FLUX" procedure, forms the 
basic double switching used to eliminate power gradients 
from the atmosphere or ground (see also Rcadhcad ct alj 
'1989| , |Myers et al. 1997| ). The ±A(j> positions are referred 
to as reference fields. 

To reduce systematic effects associated with varying 
slew times and settling of the telescope structure (particu- 
larly on the 40- meter telescope), an adjustable idle time r, 
is inserted between the A and B integrations and between 
the C and D integrations. Because the telescope does not 
move between the B and C integrations, no time delay is 
inserted there. The total duration of a FLUX procedure 
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Table 1 

Parameters for the OVRO 5.5-m and 40-m Telescopes 







5. 5- meter 


40-meter 


Center Frequency (GHz) 




31.7 


14.5 


Bandwidth (GHz) 


A// 


6 


3 


RMS Sensitivity (mK s 1 / 2 ) 


°~ rms 


1.4 


2.3 


Major Axis Beamwidth (FWHM) 


^maj 


7'. 40 ± 0'.26 a 


7'.80 ± 0'.97 a 


Minor Axis Beamwidth (FWHM) 


^min 


7'.34 ± 0'.25 a 


7'.17±0'.81 a 


Beamthrow 


A(f> 


22'. 16 


21'.50 


Main Beam Solid Angle (lCT 6 sr) 




5.21 ±0.03 


5.35 ± 0.02 


Beam Solid Angle (10~ 6 sr) 


n a 


7.92 ±0.28 


7.28 ± 0.28 


Beam Efficiency 


Vb 


0.658 ±0.024 


0.735 ± 0.029 


Aperture Efficiency 


Va 


0.476 ±0.017 


b 


Sensitivity (mK Jy — x ) 


r 


4.10 ±0.15 


21.30 ±0.82 



a Parameters are for the average of the ANT and REF beams. 

b For the under-illuminated 40-meter telescope, the physical aperture is not well determined (see tJ6l 



is thus r = 4r s + 2t;. For all of the PJNG5M observa- 
tions presented here, t s = 20 s and Tj = 10 s , so that a 
FLUX procedure typically requires 100 s . For each set of 
measurements, the quantity 



FLUX : 



i(AT B + AT c 



AT A - AT, 



D 



(2) 



is formed, with associated standard deviation (SD) esti- 
mated by summing the variances of the individual (A-D) 
integrations. 



T(<|>-A4» 



t«0 



T((j) + A(j)) 



B 
C 
D 

Fig. 2. — Double switching pattern in a FLUX procedure, shown 
here for an arbitrary temperature field on the sky. The sense of the 
beams is ANT left, REF right. 

From Figure pi it can be seen that the Dicke switching 
produces 

AT AiD =T(ci>-A4,)-T(cl>) 

AT s ,c = T(0)-T(0 + A0), (3) 

so that the double-switched FLUX is given by 

FLUX = {T(</>) - T((p + Ac/))} 

-{T(0-A0)-T(0)} 

(In Eq. |l|, Tj_|_i and Tj_i are the temperatures at the 
±A(f> positions, respectively.) If a source under observa- 
tion is smaller than the beamthrow of the telescope, i.e., 



T src O + Acp) = T SIC (<j) ~ A<j>) = 0, from Eq. | we see that 
the contribution to the FLUX is just 2T src (0), or twice the 
power increment that would be measured with a single- 
difference observation. 

On both telescopes, gain variations in the amplifiers are 
removed by referencing to noise diodes. Signals from the 
diodes are injected just behind the feed horns (see Fig- 
ure H) and are subject to the same receiver gain varia- 
tions as the astronomical signal. The diode measurement, 
known as a "CAL" procedure, is identical to the FLUX 
procedure just described, with the exception that the tele- 
scope does not move between the A and B or the C and 
D integrations. The diode remains off during the A and 
D integrations and is turned on during the B and C in- 
tegrations, so that from Eq. the recorded CAL is just 

-^diodc- 



4. RECEIVER CHARACTERISTICS 

On the 5.5-meter telescope, load switches identical to 
the Dicke switch provide an observing mode where the 
background power is input from two internal loads of tem- 
perature Tl — 20 K, providing a stable background against 
which to measure the noise diode powers even during peri- 
ods of bad weather. On the 40-meter telescope, no internal 
loads are available, and all diode measurements are per- 
formed against the sky. 

Both receivers exhibit a small degree of non-linearity, 
i.e., the measured power increment against a source in the 
presence of the typical background power level underes- 
timates the p ower increme nt at the front of the feeds by 
5 — 10% (see |Leitch 1998 for details). As a result, care 
must be exercised when comparing noise diode measure- 
ments against the internal loads to observations against 
the sky, as the background powers are in general differ- 
ent, leading to variations in the FLUX/CAL ratio as large 
as 6%. On both receivers, this effect has been measured 
and can be removed to high accuracy; at both frequencies, 
the non-linearity correction contributes <J 1% to the final 
calibration error. 
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Table 2 

Coordinates of the Ring5m Fields 



Field a(J2000) 5(J2000) 
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Table 3 

Absolute Flux Density and Temperature Scale 



Source 


31.7 GHz a 


14.5 GHz b 


Jupiter 


152 ±5 K 


175 ± 9 K 


DR 21 


20.60 ± 0.68 Jy c 


22.87 ± 1.07 Jy c 


Cas A 


164.18 ±5.45 Jy c 


313.04 ± 14.8 Jy c 


Crab 


307.15 ± 10.14 Jy c 


426.9 ±20.17 Jy c 


3C 84 




25.84 ± 1.26 Jy 


3C123 




5.85 ±0.28 Jy 


3C 218 




4.69 ±0.23 Jy 


3C 286 


2.02 ±0.07 Jy 


3.80 ±0.22 Jy 


3C 353 




8.44 ± 0.40 Jy 


NGC 7027 




6.20 ±0.30 Jy 



a Assuming a brightness temperature for Jupiter of 152 ± 5 K at 31.7 
GHz and flj up = 6.656 X 10~ 7 sr at 1 au. 

b Assuming a ratio of S DR 2i(14.5 GHz)/5 dr 21 (31.7 GHz) given by Eq. [f|. 

c Flux density as seen by the telescope beam; not an accurate measure 
of absolute flux density. 




Fig. 3. — A schematic of the 31.7 GHz receiver front end. Tal 
and Trl represent the contribution to the system temperature of 
the ANT and REF internal loads, respectively. Similarly, Tan is 
the temperature of the noise diode on the ANT side. The isolation 
of the Dicke switch (denoted D) when pointed at the ANT side is 
indicated by £a and when pointed to the REF side, by £r. 



The fast switching is accomplished by a Dicke switch — 
a three-port, wide-band circulator whose direction is de- 
termined by the polarity of the magnetic field within its 
ferrite core. With the switch directed at one input port, a 
small amount of signal from the other input port is trans- 
mitted to the output port, typically at the level of 1 — 2%, 
known as the isolation, denoted £. The directional isola- 
tions of the switch, i.e., the isolations with the switch in 
ANT or REF position, denoted £4 and £r respectively, 
depend sensitively on the impedance match at the three 
ports and in general are not equal. 

From Figure ^, it can be seen that if £4 7^ £r, even a lin- 
ear temperature gradient can partially survive the double 



switching, contributing an additional term 
AFLUX~A<^(£ R -60, 



(5) 



to Eq. 3], which, given typical values of £4 and £r, can 
be as large as 100 /iK for sky gradients as small as 
1 mK/arcminute. We believe this effect to be responsi- 
ble for year-to-year fluct uat ions in the mean observed in 
the RING5M data (see §7.4), which we expect a priori to 
be zero (see §||. 

5. FLUX DENSITY SCALE 

At 31.7 GHz, the flux density scale of the 5.5-meter tele- 
scope is based on a 31.4 GHz measurement of the bright- 
ness temperature of Jupiter, 



Tj = 152 ±5 K 



(6) 



(Dent 1972). During 1996, daily calibration of the in- 
ternal noise diodes at both frequencies was achieved by 
comparison with a set of secondary standards whose flux 
densities were measured relative to Jupiter (see Table [}]), 
using an ephemeris distance for Jupiter and assuming 
^ju P = 6.656 x 10~ 7 sr at 1 au. Ratios to Jupiter were de- 
termined at three epochs during which calibrator sources 
were observed for several days, as well as from daily obser- 
vations of the sources (see below) . Although the supernova 
remnants Cas A and the Crab nebula are partially resolved 
on both the 5.5-meter telescope and the under- illuminated 
40-meter telescope ((9 snr ~ 4'), pointing on these sources 
has proven reproducible to high precision, making them 
suitable as relative calibrators. 

Daily calibrator observations were interleaved with ob- 
servations of the RING5M fields, so that every twelfth field 
(every eight hours) was replaced by a 40-minute scan on a 
calibrator source. To avoid selective depletion of data from 
any three fields, the set of secondary calibrators was cho- 
sen so that at least one would be visible at any time, and 
calibrator scans were precessed daily by one field, result- 
ing in a uniform reduction in sensitivity of only 4% over 
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the entire Ring. At 14.5 GHz, observations of calibrator 
sources were restricted to lie within ±5° of the zenith angle 
of the Ring, and a correspondingly larger set of secondary 
calibrators was used to satisfy this condition. 

5.1. DR 21 and the 14-5 GHz Flux Density Scale 

The 31.4 GHz brightness temperature of Jupiter in Eq. || 
is based on a ht to the spectrum of the HII region DR 21, 



S v = 26.78 -5.63 log u GSz Jy, 



(7) 



with an asso ciated error of ±3% over the range 7 — 40 GHz 
(Dent 1972). Because of its location in a complicated re- 
gion of the Galactic plane, DR 21 itself is not used as an 
absolute calibrator in this experiment. Since the reference 
fields in a FLUX procedure are displaced azimuthally by 
22' (see Figure [l]), as DR 21 is tracked on the sky, emis- 
sion in the ring of radius 22' around the source rotates 
through the reference beams, making the measured flux 
density a function of the source parallactic angle ip p , de- 
fined as the anglebetween the great circle passing through 
the source and the zenith and the great circle passing 
through the source and the celestial poles (throughout 
this paper, we have folded the parallactic angle into the 
range (—90°, 90°)) (see Figure 0). Independent scans on 
DR21, however, agree over the full parallactic angle range 
to within the scatter of the data. The resultant variation 
in the apparent flux density of DR 21 can therefore be 
removed, making it suitable as a relative calibrator. 



fact can be used to establish a 14.5 GHz flux density scale 
relative to the same measurement of Jupiter at 31.4 GHz. 
The flux density scale so derived is in excellent agreement 
with the 14.5 GHz measurement of Jupite r's brightness 



temperature relative to DR 21 obtained by Gary (1974) 



using the Goldstone 64-m telescope. Therefore this scale, 
with the 5% errors reported by Gary, is the one used to 
calibrate the 14.5 GHz data. 

5.2. Cas A and Re- calibration of the Saskatoon 
Experiment 

The Saskatoon group has recently reported a measure- 
ment of intrinsi c anisotropy relative to a fit to the spec- 
trum of Cas A ( Netterfield et al. 1997| ). The power they 
detect on degree scales has led to speculation that their as- 
sumed flux density for Cas A may be too high. Since Cas 
A is used as a calibrator source in the RING5M experi- 
ment, the OVRO observations of Cas A not only provide 
a corroborative check on the Saskatoon calibration, but 
can in principle refer the Saskatoon measurements to a 
flux density scale based on Jupiter, thereby reducing their 
calibration uncertainties to a few percent. Since Cas A is 
comparable in size to the 5.5-meter telescope beam, how- 
ever, our measurement cannot be compared directly with 
the Saskatoon value, but must be multiplied by the factor 



/ = 



/ n P cas (0,0)iWM)dO' 



(8) 



34 
33 
30 
18 
16 



1.1 
1 

0.9 



m 0.8 




Fig. 4. — (Top panel) A typical scan on DR 21. The apparent 
flux density is dominated by contamination in the reference beams 
and is a strong function of parallactic angle. (Bottom) Dependence 
of the switched flux on parallactic angle, shown here in a compi- 
lation of 34 independent tracks on DR 21 from 1995. Each scan 
was separately calibrated and normalized to the flux over the ifi p 
interval (—80°, —70°). Shown also is the best fit model for the ip p 
dependence. 

It is found that the fit to Sdr 21 (Vv) derived at 31.7 
GHz also removes the parallactic angle dependence of the 
flux density at 14.5 GHz, indicating that the contaminat- 
ing flux in the reference fields has the same spectrum as 
DR 21 between 31.7 and 14.5 GHz. This is not surprising, 
since all of the emission in the region surrounding DR 21 is 
thought to belong to the same HII complex. Since the fre- 
quency dependence of the contaminating flux is separable 
from its parallactic angle dependence, the observed flux ra- 
tio at the two frequencies should follow from Eq. 0, and this 



where P cas (0,(/)) is the source brightness distribution, and 
-fovro(#, 0) is the normalized antenna power pattern, given 
by the central lobe in Figure |l|. 

The correction factor in Eq. || was determined from a 
32 GHz map of Cas A, made with the 100-m Effclsberg 
telescope (Morsi 1997, private communication), as a tem- 
plate for P cas (6,(f). The effect of uncertainties in the 5.5- 
meter telescope pointing, as well as in the determination 
of the telescope beam, was modeled via Monte Carlo simu- 
lation. Gaussian-distributed pointing positions were gen- 
erated for 10 5 events, with a a = &z = 0.'3 (typical for 
both telescopes) centered on the nominal pointing posi- 
tion; bcamwidths were drawn from Gaussian distributions 
centered on the best-fit 6 m i n and 6> ma j . The 68% confidence 
interval of the resulting distribution for / is given by 



/ = i.is: 



o.oi- 



(9) 



With the flux density as seen by the OVRO 5.5-meter tele- 
scope beam given by S' ovro = 164.18±5.45 Jy (see Table ||), 
and accounting for the secular decrease in the flux density 
of Cas A since 1994 ( Baars et al. 1977 ), application of 
Eq. I yields 



Sq 



195.5911$ Jy 



(1.05±0.04)5 sask 



(10) 



(11) 



(Netterfield 1995), consistent with the Saskatoon calibra- 
tion. 



Leitch et al. 



7 



6. CALIBRATION 

The power received from a source of specific intensity I 
uniformly filling the main beam f2 m of the telescope power 
pattern is given by 



P = ^IA p n a fl m Av, 



(12) 



where r\ a is the ratio of the effective aperture of a dish 
to its physical aperture A p , known as the aperture effi- 
ciency. The power in instrumental units is converted to 
physical units by comparison with the power emitted by a 
calibrator source of known intensity 7 ca i, typically a noise 
diode internal to the receiver. Since the radiation from an 
internal noise diode fills the beam solid angle il a of the 
telescope, the ratio P/P ca ,\ is given by 



P 



P. 



cal 



-^cal^a 



s, 



cal 



(13) 



so that the intensity for a source filling the main beam is 
given by 

,14) 

Atmospheric attenuation reduces the observed intensity of 
a source by a factor k(Z) = exp(— tA(Z)), where r is the 
atmospheric opacity at zenith and A(Z) ~ sec (Z) is the 
airmass. Thus, to recover the intensity of a source above 
the atmosphere, we must compute 



P 

Peal 



Seal 



1 



k(Z) 



(15) 



Throughout this paper, we will alternately use inten- 
sity I, brightness temperature Tb (the equivalent Rayleigh- 
Jeans (R-J) temperature of a source filling the main beam 
fl m ) and antenna temperature Ta (the equivalent R-J tem- 
perature of a source filling the beam solid angle fi„) . These 
are related simply by 



ISL 



2kT B _ 2kT Af ^ 



whence 



A 2 



rr T A 

J-B = : 

lib 



A 2 



(16) 



(17) 



where n b = Q m /Q a . 

The mean 31.7 GHz zenith atmospheric opacity at the 
OVRO during 1996 is determined by fits to the daily cali- 
brator source observations, yielding T31.7 = 0.045 ± 0.002, 
consistent with the mean annual zenith opacity estimated 
from a water vapor radiometry (WVR) system at the 
OVRO during 1994-1996. The error in 7 introduced 
by adopting a constant m ean opacity during 1994-1996 
is < 0.02% QLeitch 199f ). Extrapolation to 14.5 GHz 
of WVR opacities measured at 31.4 and 20 GHz yields 
T U .5 = 0.023. 

Maps of the main beam, shown in Figure |l|, were 
obtained in 1995 from raster-scans across Jupiter and 
3C 84 (for the 7'.4 beam, Jupiter is approximately a point 
source). At 14.5 GHz, scans were restricted to lie within 
±5° of ^Ring so that the resulting beam map is the one 
appropriate for calibration of the RING5M data (see §0). 



At both frequencies, the main beam solid angle Q m is de- 
termined to an accuracy of < 1%. 

As discussed in §|| three independent estimates of the 
noise diode flux density 5 ca i were obtained each day dur- 
ing 1996, from which we deduce that the intrinsic output 
of the diodes at 31.7 GHz varies by < 1%. At 14.5 GHz, 
instrumental effects resulted in a variation of < 2% in the 
diode power output. In the Rayleigh-Jeans regime, 



S, 



2kT, 



cal 



cal 



A 2 



a, 



(18) 



so that the flux density of the diodes can also be de- 
termined by measuring the diode antenna temperature. 
These measurements are performed in the standard man- 
ner by comparison with external loads of known temper- 
ature and form the basis of our calibration prior to 1996, 
when calibrator sources were not observed on a regular ba- 
sis. The uncertainty in the diode temperature determined 
from these measurements, when combined in quadrature 
with the error in the non-linearity correction (see §^), is 
3.2%. Note, however, that this calibration method requires 
a separate measurement of fi a from Eq. [ll| so that uncer- 
tainties in the temperature scale enter twice if Q a and T ca \ 
are not determined simultaneously. 

At 31.7 GHz, the first calibration method yields 



a?// 2 = (0.6%)^ + (3.3%)| cai , 



(19) 



for a total calibration uncertainty of 3.4% from 1996, while 
the second gives 



a] /I 2 = {0.6%)% + (3.5%)L 



(3.2%) 2 



(20) 



or a total calibration uncertainty of 4.7% prior to 1996. 
At 31.7 GHz, the mean of the calibration errors from the 
three independent seasons (see below), weighted by the 
measurement error in the variance from each season, gives 
a total calibration uncertainty of 



°"J31.7 GHz 



4.3%. 



(21) 



7. DATA SELECTION 



On the 5.5-meter telescope, a total of three seasons of 
data were obtained at 31.7 GHz from the winter of 1993 
to the spring of 1996. A typical observing season at the 
OVRO lasts from early October until mid May. On the 
40-meter telescope, construction of the 14.5 GHz receiver 
was completed in the spring of 1994, and only two sea- 
sons of data were obtained. Due to procedural differences 
in observing strategy between the first and second halves 
of 1996, however, these data are divided into two seasons 



which are analyzed separately below (for details, see Leitch 
19981). 



7.1. Miscellaneous Edits 

The first level of FLUX editing consists of rejecting all 
data taken when the receiver was saturated during peri- 
ods of high atmospheric water content. Since the receiver 
typically remains saturated during all four segments of the 
FLUX procedure (see §||), these are readily identified as 
data with standard deviations identically zero. This re- 
jects less than 1% of the data. Next, data taken during 
excessively windy conditions, leading to tracking errors, 
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are identified by the excess time taken for a FLUX pro- 
cedure to complete. Any data for which the difference 
between the actual and expected duration is > I s are ex- 
cised, eliminating between f — 10% of the data. The range 
indicates the spread between the three major divisions of 
data at each frequency, described above. 

The non-linearity correction discussed in §|4| requires in- 
terpolation of the total power onto the FLUX data. Any 
FLUX for which no bracketing power measurements are 
found within l h is rejected, typically affecting 1% or less 
of the data. No non-linearity correction is applied to the 
14.5 GHz data, as both CALs and FLUXes are measured 
against the sky, and thus against the same power back- 
ground. 

7.2. Statistical Edits 

The noise diodes are sampled every 15 minutes to re- 
move gain fluctuations from the data. The scatter in these 
diode measurements therefore provides an unbiased crite- 
rion for culling data affected by rapid gain variations on 
both telescopes, and by atmospheric fluctuations on the 
40-meter telescope. (Diodes at 14.5 GHz are measured di- 
rectly against the sky, while diodes on the 5.5-meter tele- 
scope are measured against internal loads, reducing the 
contribution of atmospheric fluctuations to scatter in the 
CALs by a factor of ~ 100.) In this step, 4 - 12% of the 
data at 31.7 GHz and 14 - 40% of the data at 14.5 GHz 
are rejected. 

Next, a series of statistical edits are applied to the 
data. These have previously been described in Myers et 
al. (1997), and the same notation is retained here for 
consistency. A combination of sliding buffer edits is im- 
plemented, where in each case we form the test statistic 



U = X i ja l 



(22) 



where the data Xi have been divided by <Ji th , the expected 
thermal noise for the i th measurement (a rms = 1.4 mK s 1 / 2 
at 31.7 GHz, and er rms = 2.3 mKs 1 / 2 at 14.5 GHz). The 
tested quantity X can be the SD or SW, a combination 
of the A-D integrations (see §p|) which cancels the signal 
in the far field, used to reject data during periods of high 



residual atmospheric or ins trumental fluctuations (see My- 



ers et al. 1997, Lcitch 1998). Successive buffers of N points 
are constructed in time, where the width of the buffer is 
constrained to be no more than two hours. For each buffer 
j, we compute the mean 



1 J+N-l 

v £ 1 



and standard deviation 



N 



j+N-l 



(23) 



N 



1/2 



(24) 



A point is rejected if there exists no buffer containing that 
point for which either tj < t max or <Xj < <r max , des- 
ignated "meanX[N, i max ]" and "sigX[N, a max ], v respec- 
tively. These filters are applied to the combined data from 
all RING5M fields, as their primary purpose is to reject 



data affected by the atmosphere, regardless of the field be- 
ing observed. In addition to these buffer edits, we employ 
a simple point-by-point filter which rejects data for which 
U > t max , designated "X[i 

max\ ■ 

A final edit, and the only edit in which the FLUX data 
themselves are used as a rejection criterion, is an itera- 
tive, field-by-field 4cr outlier rejection. This serves to re- 
ject isolated spurious signals due to local radio-frequency 
interference, and typically affects < 0.5% of the data. The 
combined edits reject 50% of the data at 31.7 GHz and 
40% at 14.5 GHz and are summarized in Table [|. 

Th e ef fect of the editing on the weighted field means 
(see |7.3| ) was investigated by reducing the data for a wide 
range of editing parameters; the mean standard deviation 
per field introduced by varying the cutoffs is found to be 
< 4 /xK (see Figure ||). 



e = 14.86 /J.K 
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Fig. 5.— Reduction of the 1994 RING5M 31.7 GHz data for 
15 editing schemes for unweighted means (top panel), and the same 
reductions using weighted means (bottom panel). Even in the case 
of the unweighted editing schemes, the mean standard deviation per 
field (shown at the top right of each plot) is smaller than the statis- 
tical error bars, and in the case of the weighted data, is practically 
negligible. 

At 31.7 GHz, measurements of the noise diodes against 
the internal loads can be used to estimate the residual con- 
tribution of the atmosphere to the filtered data standard 
deviations; while the distribution of diode SDs is consis- 
tent with the expected thermal noise, the distribution of 
FLUX SDs against the sky peaks at 1.4 times the ther- 
mal limit, with a considerable skew to higher SDs. We 
attribute this excess noise to residual atmospheric fluctua- 
tions not removed by the fast switching (recall that on the 
5.5-meter telescope, the ANT and REF beams depart sig- 
nificantly beyond ~ 0.5 km, while the typical scale height 
of water vapor is ~ 2 km). At 14.5 GHz, the peak of the 
FLUX SD distribution is within 10% of the thermal limit, 
consistent with the lower atmospheric opacity at 14.5 GHz 
and the better overlap of the 40-meter telescope beams. 

7.3. Field Means 

Data for each RING5M field are acquired in 40-minute 
scans, during which ~ 21 double-switched FL UXe s are col- 
lected. After the data filtering described in §7.1 -§7.2, the 
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Table 4 
Data Edits 



31.7 GHz 14.5 GHz 



Edit a 


1994 


1995 


1996 


1995 


1996.1 


1996.2 


Saturated (%) 


0.05 


0.77 


0.15 


0.20 


0.64 


0.03 


Outlier (%) 


0.46 


0.24 


0.14 


0.21 


0.15 


0.37 


No Powers (%) 


0.69 


1.41 


0.00 


0.00 


0.00 


0.00 


SD[2.5] (%) 


1.87 


1.77 


2.34 


2.34 


1.92 


2.72 


sigSW[25,7.5] (%). . 


2.25 


1.67 


2.07 


0.02 


0.00 


0.00 


CALs (%) 


3.99 


4.06 


11.39 


24.61 


39.80 


13.58 


Excess time (%) . . . 


9.72 


1.00 


3.92 


24.60 


1.29 


3.67 


meanSD[25,2.0] (%) 


22.64 


36.87 


35.53 


3.16 


1.11 


1.38 


Total rejected (%) . 


41.67 


47.78 


55.53 


55.13 


44.91 


21.75 


Total 1 ? 


68,359 


98,617 


107,494 


32,743 


29,209 


40,584 


Note. — Percentage fi 


gures are 


percent of data points that 


were rejected. 







a For an explanation of notation, see §|7j. 
b Total number of data points before editing. 



weighted mean A7\ ± for each field i is computed as 
follows: 



1 



2 _2 



w. 



Ni 
Ni 

i=i 



3=1 



(25) 
(26) 
(27) 

(28) 

(29) 



where ATy is the double-switched temperature from a sin- 
gle FLUX procedure (see §0), and Ni is the total num- 
ber of FLUXes recorded for field i. Formally, the weights 
Wy should be chosen so that the sample mean is a maxi- 
mum likelihood estimator for the mean of the underlying 
distribution, i.e., Wy = cr^ 2 - In general, however, the 
standard deviation reported with each datum (typi- 
cally 1.4 x the thermal noise at 31.7 GHz) underestimates 
the scan standard deviation by a factor of 2-3, presum- 
ably due to atmospheric fluctuations on timescales longer 
than a single FLUX measurement. The scan variance <r|. 
is thus a better estimate of the real error, and we take 
w 



c 2 + erf,- __) 1 when computing statistics for the 
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Fig. 6.— PJNG5M results for 1994-1996. (Top panel) Three 
years of data, with no means subtracted. Subtraction of the 1994 
and 1995 data sets (middle panel) shows that the Ring has changed 
by a constant offset (AT = 119.41 ^K, \l = 1-34), while the 1996 
and 1995 data sets are consistent with no change (AT = —9.93 /^K, 
Xr = 0.71) (bottom panel). 



7.4. Mean Levels 

The 31.7 GHz RING5M field means for 1994 - 1996 are 
shown in Figure ^. Between 1994 and 1995, the mean 
level of the Ring changed by approximately 120 /iK, while 
between 1995 and 1996, the difference is consistent with 
zero. Subtraction of the data sets shows that in each case 
the shift is consistent with an offset which is constant from 
field to field and thus does not affect our estimate of the 
sky variance (a variable point source contributes to the 
measured signal in the first field (see §[lO]), which affects 
the two neighboring fields through the double switching). 
The constancy of these offsets implies that they are instru- 
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mental in origin. 
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Fig. 7.— Mean levels for the 31.7 GHz Ring data, 1994-1996, 
where bins contain equal numbers of data points. The large step 
in the mean level coincides with a HEMT change on day 91, 1994 
(indicated byJJie dotted line). Assuming variations in the mean are 
given by Eq. Q the predicted mean of the Ring at the two epochs 
in 1994 and 1995 when measurements of the directional isolations 
were made are indicated by open triangles, normalized to the mean 
observed in 1993. 

Figure |?] shows that the largest shift in the mean level 
occurred after a HEMT changeover on day 91 of 1994. 
In §||, we showed that a mismatch in the directional iso- 
lations of the Dicke switch can lead to incomplete can- 
cellation of ground or atmospheric temperature gradients. 
Replacement of the HEMT will undoubtedly change the 
impedance at the output port of the switch, quite likely 
resulting in a change in isolation. At 31.7 GHz, three 
measurements were made of the isolations, one before the 
HEMT changeover, and two after. We find that when the 
error term given by Eq. || is normalized to the RING5M 
mean in 1993, the mean levels predicted by the isolations 
measured in 1994 and 1995 (shown as open triangles in 
Figure ^) are in excellent agreement with the data; in 
each case the predicted and measured mean levels agree 
to within errors and imply a constant temperature gra- 
dient of A0^- = — 16 mK across the beamthrow of the 
telescope. 

8. CONCORDANCE 

Edited and calibrated means for the 36 fields, represent- 
ing a combined total of approximately 4,500 hours of data, 
are shown in Figure |l A mean level has been subtracted 
from each season. Within each season, we test for inter- 
nal consistency by dividing the data set into two halves 
in time. While the time spanned by a given season is in- 
sufficient to permit a sensible comparison of data taken at 
night to data taken during the day (see §[?]), the data are 
nevertheless quite repeatable on timescales over which the 
field positions have precessed significantly relative to the 
Sun; linear correlations (Pearson's r) between halves fall 
in the range r bs = 0.72 — 0.89. 

At both frequencies, the data show structure well above 
the noise and repeatable from year to year. At 31.7 
GHz, the season-to-season deviations in the field means 
are Gaussian, and correlations between seasons are r b s > 
0.89. The probability of observing correlations this high 
under the hypothesis that the sky temperatures are uncor- 



related is p(r > r b s ) < 6 x 10~ 8 . Under the hypothesis 
that the data are completely correlated, as we should ex- 
pect if the signals are dominated by the microwave back- 
ground, these correlations fall within the 68% confidence 
region for r b S) given typical season-to-season field errors 
of ~ 20 /iK (antenna temperature) . Because of the double 
switching, the RING5M field means are not statistically 
independent; all probabilities quoted in this paper take 
the effect of the switching into account. In addition, we 
assume that the unswitched sky temperatures are drawn 
from a Gaussian distribution. 
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Fig. 8.— RING5M 31.7 GHz (top) and 14.5 GHz (bottom) field 
means (antenna temperature) from winter of 1993 (7^.94) to winter 
of 1995 (lige). The rms is denoted a, and e is the mean error per 
field. At bottom left is the beam pattern produced by the double 
switching, i.e., the effective point-spread function for the OVRO 
telescopes. 

At 14.5 GHz, the agreement is also good, with the ex- 
ception of the RING5M field at 21 h 04 m (OV5M2104), in 
which an anomalously high (and to date unexplained) sig- 
nal was seen during the first half of the 1995 season, and 
OV5M0024, which is dominated by a bright variable point 
source (see §[l0|). The two fields adjacent to OV5M0024 are 
affected by the same point source due to the double switch- 
ing. These effects collectively reduce season-to-season cor- 
relations at 14.5 GHz to the range r D b s = 0.72 — 0.87. 

9. TERRESTRIAL CONTAMINATION 

Since the features we detect at both frequencies are fixed 
in sidereal time, local radio frequency interference (RFI) 
is an unlikely explanation for the structure observed in 
the RING5M data. Nonetheless, if there is ground-based 
interference, the change in telescope elevation as a field 
is tracked through transit will introduce a characteristic 
parallactic an gle dependence into the data for any one 
field (see §5.1), and we can look for this signature in the 
RING5M data. 

The means for all 36 fields are shown in Figure ^[ binned 
in parallactic angle. We find no parallactic angle depen- 
dence in the 31.7 GHz data, indicating that these data are 
free from RFI contamination. The 14.5 GHz data, on the 
other hand, show a large variation in amplitude with par- 
allactic angle, indicating significant contamination. This 
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pattern, however, occurs with the same amplitude in each 
field (Figure ^ is a compilation of data from all 36 RING5M 
fields), demonstrating that RFI contributes only a mean 
level, constant over long timescales, so that our measure- 
ment of the variance should be unaffected. 
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Fig. 9.— Co-added 31.7 GHz (top) and 14.5 GHz (bottom) data 
from all 36 RING5M fields, binned in parallactic angle. The 14.5 
GHz data show the characteristic signature of interference from the 
ground, while the 31.7 GHz data appear free of contamination. 

Although the parallactic angle dependence is a good in- 
dicator of the robustness of our data to interference, we can 
definitively rule out terrestrial contamination as a source of 
structure at 14.5 GHz by observing the fields at lower cul- 
mination. Since this changes the position of the telescope 
beams only relative to the ground, signals from the sky will 
remain unaltered, while any time-dependent fluctuations 
from the ground will be shifted by 12 hours. This test 
was performed on 1996 Sep 16-30. (The strength of the 
signals at 14.5 GHz permits a reasonable detection in this 
relatively short time.) The comparison of the mean upper 
culmination 14.5 GHz data set with the lower culmination 
data is shown in Figure |l^. Field errors are estimated by 
reducing the 1996 upper culmination data in 2- week sub- 
sets; the error in the mean for each field is typically about 
62 fiK in antenna temperature. 



beam pattern, which is inverted at lower culmination rel- 
ative to the upper culmination beam. 

10. DISCRETE SOURCES 

Mosaicked VLA observations covering the RING5M 
fields to the 3% contour of the 7'.4 beam were made in 
1994, with an rms sensitivity at 8.5 GHz of 0.21 mJy. The 
sensitivity of the 5.5-meter telescope to point sources is 
4.1 mK Jy _1 , so that these observations allow detection 
of any source contributing J> 13 (4tr) to our highest fre- 
quency data with an accuracy of <; 3/zK (la), assuming 
a < +1 (where S v oc v a ). During the period bracket- 
ing the 1996 RING5M season, we obtained multi-epoch 
VLA observations at 8.5 and 15 GHz of the 39 sources (56 
discrete components in all) found in the original 8.5 GHz 
survey. These observations allow extrapolation of source 
flux densities to 31.7 GHz and removal of variable contri- 
butions to the RING5M data from sources which vary on 
timescales 1 month. 

From October 1995 - May 1996, RING5M sources were 
observed at 11 separate epochs, in BnA, B, CnB, C and 
DnC configurations. Due to the time constraints of the 15 
GHz observations (the approximate 10 minute sensitivity 
of the VLA at 15 GHz is 0.17 mJy, compared to 0.045 mJy 
at 8.5 GHz, so that 80% of our observing time was spent at 
15 GHz), not all sources could be observed at each epoch. 
Sources which early-on showed significant variability were 
observed at each epoch, while any remaining sources which 
could not be observed at a given epoch were observed dur- 
ing the next. Typical rms noise in the maps was 0.25 mJy 
at 8.5 GHz and 0.5 mJy at 15 GHz. 
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Fig. 10. — 14.5 GHz upper culmination field means (dot-dashed 
line), shown with the lower culmination data (solid line). The good 
agreement of these two data sets effectively rules out terrestrial RFI 
as the source of the observed structure in the Ring. The large differ- 
ences in the first, second and last fields are due to discrete sources 
(see text). At bottom left is the effective beam pattern for double- 
switched observations. 

The lower culmination data show the same structure, to 
within errors, as the upper culmination data, demonstrat- 
ing that the origin of these signals is undoubtedly celestial. 
The large differences in fields OV5M0024, OV5M0104 
and OV5M2344 are due in part to variability of a point 
source which dominates the signal in field OV5M0024, 
and to the slightly asymmetric 14.5 GHz double-switched 
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a 

Fig. 11. — The histogram of spectral indices for sources selected at 
8.5 GHz, derived from the 47 discrete source components detected 
at both 8.5 and 15 GHz. Here, p(ct) is the number of sources with 
spectral index a ± 0.15. 

For each source, the visibilities from all epochs at which 
the source was observed are combined to form a single 
high-sensitivity map at each frequency from which an ac- 
curate model for the source spatial structure can be deter- 
mined. Several sources have multiple components, consist- 
ing of a core and one or two prominent lobes, while oth- 
ers show two point-like components closer than one would 
expect from random superpositions of unrelated sources. 
Model-fits to these combined maps also provide a measure- 
ment of the mean flux density at each frequency. These 
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are given m Table g for sources detected at both frequen- 
cies, along with the source flux densities extrapolated to 
31.7 GHz. 

Although sources detected at 8.5 GHz but not at 15 GHz 
- nine in all — have falling spectra and will contribute 
negligibly to the RING5M data, we can take a Bayesian 
approach (Sivia 1996) to estimate their flux densities at 
higher frequencies. The distribution of spectral indices for 
sources selected at 8.5 GHz, shown in Figure [ll], is con- 
structed from the 47 components detected at both 8.5 and 
15 GHz. Once the distribution of spectral indices is known, 
the probability distribution of the 15 and 31.7 GHz flux 
densities, given the measured flux density at 8.5 GHz and 
an upper limit at 15 GHz, can be constructed for each 
source via Monte Carlo simulation, yielding the maximum 
likelihood estimates and corresponding 68% confidence in- 
tervals given in Table |. 

The estimated contribution of point sources to the 1996 
31.7 and 14.5 GHz data is shown in Figure [l^. Few of 
the sources are bright enough or close enough to a field 
center to contribute a significant signal, with the notable 
exception of a ~ 100 mJy source which dominates the sig- 
nal in field OV5M0024 and affects the two adjacent fields 
through the double switching. If we exclude OV5M0024 
and its flanking fields, the rms at 31.7 GHz due to discrete 
sources is 12 /iK (antenna temperature), or < 4% of the 
observed variance, and 58 fiK at 14.5 GHz, or < 10% of 
the variance observed there. 



&A = <*z = 0'3. At 31.7 GHz, the largest uncertainty due 
to pointing is 17 /iK, while at 14.5 GHz, the largest error is 
70 /iK, comparable to the variance of the brightest variable 
source in the Ring, making pointing the dominant uncer- 
tainty. We assume this to be true even for years when the 
sources were not monitored. 

In all subsequent analysis, a single combined data set is 
used at each frequency. Although the errors due to point- 
source subtraction are not comparable from year to year, 
we can define a mean error for each field 



where Uj is the error from source subtraction from sea- 
son j, and uij is the sum of the weights from season j, 
i.e., 

w i = E ( 31 ) 

k 

where tjk is the error in an individual FLUX measure- 
ment. To these we add in quadrature the statistical errors 
for each field mean and the error per field due to pointing 
uncertainties. 

The mean source-subtracted R-J temperatures ATj and 
associated uncertainties e, at each frequency are given in 
Table § 
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Fig. 12. — Point source contribution to the 31.7 GHz Ring data 
(top), and to the 14.5 GHz data (bottom), extrapolated from VLA 
monitoring at 8.5 and 15 GHz. Error bars reflect uncertainties in 
the source fluxes due to pointing only. At bottom left is the effective 
response to a point source near a Ring field center. 

Because sources were monitored during 1996, formal 
errors from point source subtraction for this season are 
typically a few /iK. Since many of the RING5M sources 
lie on the exponential cutoff of the 7'.4 beams, however, 
small errors in pointing can produce large changes in the 
antenna temperature produced by a source, so that the 
variance due to pointing errors contributes significantly 
to the uncertainties. Field errors in Figure [j^ are the 
computed standard deviations due to pointing, assuming 
Gaussian azimuth and zenith angle pointing errors with 
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Fig. 13. — Mean source-subtracted 31.7 GHz (solid line) and 14.5 
GHz_£dotted line) RING5M data, to equal intensity scales. See 
Eq. hq for the conversion from AI (kjy/sr) to ATb (K). 



11. FOREGROUNDS 

Mean (1993-1996) source-subtracted data sets at each 
frequency are shown in Figure [l3|. The correlation be- 
tween frequencies is remarkably high, with r D b s = 0.81 and 
probability of observing a higher correlation than this un- 
der the hypothesis that the data are uncorrelated given by 
p(r > r bs) = 2.7 x 10~ 6 . The good agreement between 
several years of data from two independent instruments 
demonstrates that the observed structure in the RING5M 
data must originate outside the telescopes. Furthermore, 
as discussed in §|[ observations of the RING5M fields at 
lower culmination rule out interference from the ground. 
We are therefore confident that the signals are celestial in 
origin. 
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Table 5 

Summary of Source Flux Densities 



a (J2000) 


S (J2000) 


Ss.5 GHz (nijy) 


Si ^ c^vfr, (mjvl a 


S31.7 ghz ( m Jy) a 


00:22:12.97 


87:53:56.64 


2.10 ± 0.10 




0.20l°;^ 


00:32:41.48 


87:50:43.61 


89.86 ± 0.56 b 


107.08 ± 0.49 


126.44 ± 0.60 


02:22:09.58 


87:52:18.73 


4.14 ± 0.10 


3.67 ± 0.26 


3.12 ± 0.27 


02:31:45.60 


87:49:10.76 


4.49 ± 0.10 


5.03 ± 0.34 


5.85 ± 0.47 


04:28:27.47 


87:58:32.65 


17.89 ± 0.15 


9.75 ± 0.31 


4.33 ± 0.15 


06:30:29.55 
06:30:29.55 


87:53:18.41 
87:53:19.64 


1.84 ± 0.11 
0.66 ± 0.11 


0.69+"^ 

— 0.06 

0.32+"?!; 

— 0.10 


0.18+"" 

w u — 0.04 

0.09 + „-^ 

— 0.05 


06:34:14.18 


87:56:23.83 


2.78 ± 0.11 


1 01+ 014 

■"■■"-0.06 


°- 26 -0.04 


07:02:43.09 


87:45:09.00 


7.67 ± 0.12 b 


7.84 ± 0.40 


8.07 ± 0.46 


08:56:31.32 


87:47:58.19 


1.70 ± 0.11 


2.10 ± 0.28 


2.79 ± 0.53 


09:38:58.18 


87:50:52.59 


2.23 ± 0.12 


U - 84 -0.09 


22 + U IU 


10:18:52.68 
10:20:12.52 
10:21:12.46 


87:56:08.52 
87:56:11.76 
87:56:00.62 


8.32 ± 0.09 
5.51 ± 0.09 
22.67 ± 0.09 


5.97 ± 0.42 
4.87 ± 0.37 
17.13 ± 0.41 


3.83 ± 0.27 
4.13 ± 0.37 
11.78 ± 0.27 


10:56:50.10 


87:47:46.16 


5.19 ± 0.10 


5.97 ± 0.33 


7.76 ± 0.48 


11:46:15.05 
11:46:13.76 


87:54:55.13 
87:54:56.52 


4.88 ± 0.08 
0.99 ± 0.08 


3.59 ± 0.29 
0.39+"-" 

— . 5 


2.38 ± 0.24 


11:48:37.64 
11:48:39.63 


87:42:05.72 
87:42:09.74 


13.13 ± 0.12 b 
2.62 ± 0.12 b 


5.86 ± 0.30 
0.87 ± 0.30 


1.99 ± 0.12 
0.20 ± 0.12 


11:53:24.39 


87:56:06.25 


7.07 ± 0.10 


3.94 ± 0.30 


1.80 ± 0.17 


12:11:50.52 


87:50:54.44 


11.88 ± 0.32 


12.72 ± 0.31 


13.94 ± 0.39 


12:16:17.38 


87:51:24.27 


20.45 ± 0.23 b 


21.36 ± 0.31 


18.76 ± 0.31 


12:55:57.99 


87:48:00.82 


3.72 ± 0.10 


3.66 ± 0.27 


3.58 ± 0.33 


13:00:16.65 


87:45:09.98 


9.16 ± 0.11 


5.10 ± 0.30 


2.33 ± 0.16 


13:41:49.23 


87:48:20.03 


1.89 ± 0.09 


2.99 ± 0.29 


5.52 ± 0.70 


14:27:16.54 


87:47:39.89 


4.70 ± 0.16 b 


3.18 ± 0.26 


1.89 ± 0.20 


14:33:37.94 


87:51:07.58 


2.91 ± 0.11 




°- 27 -0.04 


15:00:10.90 
15:00:10.12 
14:59:50.43 


87:50:53.45 
87:50:57.62 
87:50:07.46 


6.17 ± 0.13 
3.91 ± 0.13 
4.68 ± 0.13 


3.74 ± 0.61 
3.30 ± 0.61 
4.21 ± 0.72 


1.92 ± 0.38 
2.63 ± 0.60 
3.65 ± 0.63 


15:02:54.05 


87:58:44.39 


3.07 ± 0.12 


1.90 ± 0.38 


1.00 ± 0.28 


15:11:29.40 
15:11:34.96 


87:55:43.30 
87:55:46.31 


1.91 ± 0.11 
1.24 ± 0.11 


2.29 ± 0.39 
n 4Q+ ' 21 

— 06 


2.92 ± 0.68 

14+ 017 
— 0.05 


15:42:38.53 


87:55:38.92 


4.39 ± 0.10 


3.30 ± 0.31 


2.25 ± 0.26 


17:03:23.72 


87:45:10.68 


5.01 ± 0.09 


4.61 ± 0.30 


4.12 ± 0.32 


17:30:38.61 
17:31:05.73 
17:30:37.89 


87:54:12.72 
87:54:12.77 
87:54:16.49 


6.12 ± 0.13 
2.46 ± 0.13 
4.85 ± 0.13 


2.21 ± 0.25 
1.37 ±0.25 
3.79 ± 0.26 


0.57 ± 0.09 
0.63 ± 0.18 
2.73 ± 0.22 


17:31:04.27 


87:56:33.13 


2.19 ± 0.24 


n fiQ+U.M 
U.83_ Q 10 


23 +U ^ 1 


17:39:48.20 


87:49:52.79 


4.05 ± 0.08 


3.99 ± 0.31 


3.91 ± 0.37 


17:45:23.55 


87:44:33.37 


15.29 ± 0.13 b 


11.26 ± 0.28 


7.48 ± 0.20 


19:00:58.12 


88:01:38.76 


6.80 ± 0.12 b 


3.59 ± 0.36 


1.53 ± 0.19 


19:29:34.27 


87:55:02.19 


3.07 ± 0.09 


3.96 ± 0.27 


5.56 ± 0.47 


19:41:43.37 


87:46:35.91 


2.55 ± 0.11 


3.21 ± 0.38 


4.37 ± 0.66 


20:32:28.60 
20:32:49.85 
20:32:52.51 
20:32:30.54 


87:59:24.91 
87:59:31.15 
87:59:25.76 
87:59:30.94 


13.96 ± 0.15 
3.79 ± 0.15 
2.87 ± 0.15 
2.38 ± 0.15 


6.57 ± 0.28 
2.24 ± 0.28 
1.19 ± 0.28 
2.01 ± 0.28 


2.40 ± 0.12 
1.11 ± 0.19 
0.37 ±0.14 
1.60 ± 0.32 


21:41:01.80 
21:41:53.75 
21:41:52.32 
21:41:03.45 
21:41:41.22 


87:58:10.62 
87:57:51.41 
87:57:53.02 
87:58:10.02 
87:57:55.77 


3.50 ± 0.07 
1.65 ± 0.07 
4.31 ± 0.07 
5.59 ± 0.07 
2.02 ± 0.07 


1.45 ± 0.27 
1.54 ±0.30 
1.98 ± 0.29 
3.80 ± 0.27 
0.99 ± 0.28 


0.45 ± 0.12 
1.40 ± 0.35 
0.70 ± 0.13 
2.27 ±0.19 
0.38 ± 0.17 


22:28:10.05 


87:50:37.98 


9.02 ± 0.09 


8.82 ± 0.30 


7.98 ± 0.31 


23:16:20.13 


87:49:41.46 


7.71 ± 0.99 


3.97 ± 0.41 


1.63 ± 0.25 


23:17:25.40 


87:52:58.42 


17.73 ± 0.15 b 


13.11 ± 0.30 


8.77 ± 0.22 



a Asymmetric errors indicate a source n °t detected at 15 GHz and is 
a maximum likelihood estimate (sec §|l0|) . 

b Variablc source. 
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Table 6 

Source-Subtracted Ring Field Means 3 - 



r iciu 


i-* 1 31.7 (h 1 ^) 


e 31.7 IM^J 


AT, , r C uK\ 
<-±-L 14.5 \h LIX ) 


£14.5 {p- 


K'V 

1VJ 


OV5M0024 


-67.78 


36.88 


-710.82 


101 


.57 


OV5M0104 


27.95 


32.78 


397.69 


96 


.38 


OV5M0144 


-33.59 


17.15 


-94.52 


20 


.86 


OV5M0224 


90.73 


15.94 


339.78 


21 


.75 


OV5M0304 


-75.61 


16.40 


-168.44 


17 


.31 


OV5M0344 


-21.54 


17.35 


-23.06 


18 


.68 


OV5M0424 


79.59 


16.14 


-78.56 


24 


.04 


OV5M0504 


-15.49 


16.32 


-56.43 


18 


.36 


OV5M0544 


-13.89 


16.51 


47.66 


16 


.77 


OV5M0624 


95.96 


15.75 


75.63 


18 


.99 


OV5M0704 


126.83 


15.91 


-89.17 


17 


.21 


OV5M0744 


13.87 


15.31 


-7.65 


17 


.91 


OV5M0824 


-124.12 


13.88 


-357.49 


15 


.46 


OV5M0904 


140.53 


14.57 


329.81 


15 


.61 


OV5M0944 


46.24 


15.79 


-174.36 


16 


.63 


OV5M1024 


19.74 


15.01 


86.55 


32 


.74 


OV5M1104 


-7.82 


15.29 


262.51 


19 


.28 


OV5M1144 


-56.77 


15.92 


-381.74 


27 


.99 


OV5M1224 


51.82 


17.35 


340.51 


45 


.91 


OV5M1304 


-40.51 


15.47 


-50.77 


25 


.38 


OV5M1344 


-102.35 


15.40 


-208.90 


18 


.50 


OV5M1424 


-8.84 


15.58 


107.09 


26 


.39 


OV5M1504 


58.72 


16.95 


225.55 


35 


.13 


OV5M1544 


-34.40 


16.27 


-65.48 


19 


.18 


OV5M1624 


-67.51 


15.86 


-51.07 


18 


.61 


OV5M1704 


-38.27 


15.76 


-61.29 


18 


.11 


OV5M1744 


13.11 


16.33 


201.98 


20 


.53 


OV5M1824 


-114.55 


15.50 


-170.47 


19 


.79 


OV5M1904 


68.41 


16.03 


105.79 


17 


.35 


OV5M1944 


-48.60 


15.31 


-121.67 


19 


.37 


OV5M2024 


-191.67 


16.51 


-297.43 


19 


.70 


OV5M2104 


176.57 


16.17 


266.81 


20 


.34 


OV5M2144 


-21.59 


17.04 


152.51 


32 


.31 


OV5M2224 


163.50 


16.90 


231.69 


24 


.88 


OV5M2304 


-134.58 


17.65 


-207.03 


18 


.50 


OV5M2344 


45.92 


27.35 


204.82 


47 


.56 



a Equivalcnt R-J temperature. 

b Errors are la rms in the sample mean. 
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Since the CMBR specific intensity is approximately 
^cmb — 2fcT cm b^ 2 /c 2 , signals from the microwave back- 
ground should be reduced at 14.5 GHz by a factor of ~ 5 
relative to those at 31.7 GHz. As can be seen in Figure [l3| , 
however, many of the RING5M fields, notably in the re- 
gions h — 3 h and 12 h — 18 h , have equal intensities at the 
two frequencies, suggesting that the emission mechanism 
may be thermal bremsstrahlung. On the other hand, the 
regions 4 h — 8 h and 20 h — 23 h show the spectral signature 
of the CMBR. 

The nature of the steep-spectrum (in temperature) sig- 
nals seen in the RI NG5M data has been investigated by 
Leitch et al. (1997) (hereafter LI). There we used a maxi- 
mum likelihood test (described in detail in §[l2|) to rule out 
contamination by flat-spectrum foregrounds under the as- 
sumption that the CMBR is observed in the presence of a 
single foreground. We model the RING5M field means in 
the Rayleigh-Jeans regime as 



ATj obs — ATj cmb -I- ATi 



i fore ) 



(32) 



where ATif ore oc v$ . Given two frequencies v\ and v 2l we 
can eliminate the ATi foro and solve for the CMBR com- 
ponent in each field as a function of the unknown spectral 
index 0; 



AT, 



icmb 



03) 



AT iobB (i/i)i/i 







AT obs (jA,)i/ 2 



a{y 1 )v 1 - a{v 2 )v 2 



(33) 



where a(v) is a factor which corrects for the R-J ap- 
proximation to a true blac kbody spectrum (at 31.7 GHz, 
a = 0.974; see Leitch 1998). As can be seen in Figure O , 
foregrounds with temperature spectral indices (3 > —1.7 
can be ruled out at the 3cr level. 



the WENSS survey ( Rengelink et al. 1997 ) show no de- 
tectable signals (see Figure |l5|). As discussed in LI, we can 
use the rms in the 325 MHz map to rule out contamination 
of the RING5M data by any foreground with temperature 
spectral index (3 < —2.2. 
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Fig. 15. — 14.5 GHz (dot-dashed line) and convolved 325 MHz 
map of the NCP (solid line) from the WENSS survey, extrapolated 
to 14.5 GHz assuming /3 = —2.7. Standard deviations are quoted for 
overlap region only (WENSS data are missing for fields OV5M1224 
- OV5M1744). The comparison demonstrates that the structure 
observed at 14.5 GHz cannot be due to steep spectrum synchrotron 
emission. 
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Fig. 14.— Likelihood function £(AT cmb ./3) (see §|l2j) for the 
RING5M data, assuming CMBR + single foreground oi tempera- 
ture spectral index f3. Plotted are the 68%, 95% and 99% highest 
probability density (HPD) intervals. 

Although the RING5M data alone cannot provide much 
discrimination among the steepest spectral indices in a 
two-component model, since the CMBR component in- 
creasingly dominates for steeper foreground spectral in- 
dices, we can use low-frequency maps of the NCP to de- 
t ermi ne limits on the observed contribution of synchrotron 
emission to the RING5M data. At 325 MHz, maps from 



If the 14.5 GHz structure is due to thermal brems- 
strahlung emission, we would expect a considerable Ha 
signature, unless the temperature of the emitting gas is 
T e S> 10 4 K; lack of any detectable structure in Ha maps 
of the NCP in fact restricts the temperature of this com- 
ponent to > 10 6 K (see LI for details). In light of this, 
flat-spectrum synchrotron or non-thermal emission from 
dust (see below) may be a more likely explanation for these 
anomalous signals. 

In LI, we also reported a significant correlation between 
the steep-spectrum foreground at 14.5 GHz and IRAS 
100 /im emission near the NCP — independent confirma- 
tion that the structure observed in the Ring is celestial in 
origin (see Figure |l6|). Drainc and Lazarian (1998) have 
recently suggested that non-thermal emission from spin- 
ning dust grains could produce emission with an apparent 
free-free spectrum, while naturally accounting for the ob- 
served correlation. 

Thermal emission from the dust itself is not expected 
to be a serious contaminant at these frequencies (fits to 
the COBE DMR data on 7° scales suggest that Od us t < 
2 /iK at 31.4 GHz, while the DIRBE spatial template in- 
dicates that the power spectrum of the dust is fal ling as 
V(V) oc l~ 3 with decreasing angular scale (I oc O^ 1 ); Kogut 



et al. 1996) 
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Table 7 

Rms for Mean 1994-1996 31.7 and 14.5 GHz Data 





C31.7 GHz a 


C14.5 GHz a 


Raw 


94.65 fiK 


249.93 


After point source subtraction 


85.41 ,iiK 


238.06 nK 


After foreground fitting, with /3f or c = — 2.2 


82.51 fiK h 


234.37 ^K c 



a R-J Temperature unless otherwise noted. 
b CMBR Brightness Temperature. 
c Extracted foreground component. 



« b a = 82.51 ijK e = 22.88 /J.K 
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Fig. 16. — The extracted CMBR physical brightness temperatures 
(at 31.7 GHz) (top panel), and the extracted foreground equivalent 
brightness temperature (at 14.5 GHz) (bottom panel), assuming a 
spectral index /3 = —2.2. Note that at 14.5 GHz, even the regions 
T - 8 h and 20 h - 231, in which the CMBR dominates the combined 
signal (see Figure h.3|), now show a noticeable correlation with the 
IRAS 100 fim intensities. 



a scaled IRAS template might be a viable method of ex- 
tracting the CMBR component, we feel that this is not 
justified; CMBR plus a single foreground with the IRAS 
spatial template is a bad fit to our multifrequency data 
(best fit xf — 10), as can also be seen by the large residual 
differences in the lower panel of Figure |l6|. Furthermore, 
other methods which lead to more restrictive limits on the 
CMBR amplitude make unjustifiable assumptions about 
the distribution of the foreground. We therefore take the 
extraction outlined in the previous section to be the most 
conservative, as it makes only the assumption that the fre- 
quency depedence of the foreground can be modeled by a 
simple power law. 

The approach we take to estimating the CMBR variance 
is the standard approach o f Maximum Likelihood estima- 
tion (see, e.g., Sivia 1996). We assume that the 36 field 
means are drawn from a Gaussian distribution 3 on the sky, 
so that the joint probability of the data set is given by 



C = 



exp(-it T C- 1 t) 
(2tt) n / 2 |C| 1/2 



(34) 



where t is the data vector, with elements AT 1 ; given by 
Eq. |3^, and C is the associated 36 x 36 covariance matrix. 
In general, the elements of C can be written 



C 



% 3 cmb 5 



(35) 



where the <r|- are the temporal covariances from the data 
(the diagonal elements <j\ are just the variances of the 
field means). The C, 



V cmb 



12. INTRINSIC ANISOTROPY 

In the foregoing discussion, we demonstrated that the 
signals detected in the RING5M experiment are consis- 
tent with a combination of steep spectrum (/3 ~ —2.2) 
and blackbody emission (/3 ~ 0), the former contributing 
97% of the variance at 14.5 GHz, the latter responsible for 
88% of the variance at 31.7 GHz. Low-frequency maps of 
the NCP were used to rule out contamination by (3 < —2.2 
emission, leaving us with a foreground which cither is free- 
free, or has very nearly the same spectral dependence, in 
either case justifying the assumption of a single foreground 
and allowing separate reconstruction of the CMBR and 
foreground components from linear combinations of the 
14.5 and 31.7 GHz data (see Figure |l|). 

Although the unexpected correlation of the 14.5 GHz 
data with Galactic IR cirrus suggests that subtraction of 

3 Note that the Gaussian approximation is equivalent to retaining the first three terms in the Taylor expansion of any arbitrary distribution 
about its maximum. 



describe the predicted spatial 
variance from the CMBR, and in general are given by the 
2-pt correlation of the effective antenna pattern convolved 
with the sky temperature field (see §|l3|). 

Once the likelihood function £(AT cm b, (3) is constructed 
for the CMBR component, our estimate of the CMBR 
variance AX" c 2 mb = C*M Cm b is obtained by maximizing C 
with respect to AT cm b. Implicit in this construction is 
the assumption of a spectral index (3 for the anomalous 
foreground; note however, that even if all of the signal at 
14.5 GHz were foreground, the relative contribution of this 
foreground to the total variance at 31.7 GHz is less than 
25% for spectral indices as flat as (3 — —2.2 (see Table 0). 
Thus our estimate of the CMBR variance is not strongly 
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dependent on the exact choice of foreground spectral in- 



dex, as can also be seen from Figure 14 



A conservative approach is simply to integrate out the 
dependence of the likelihood on (3, i.e., 



£(AT cmb ) = J £(AT cmb ,p)d8, 



(36) 



where we restrict the integration to bounds defined by a 

reasonable prior e g one uniform for — 3 < ft < 9. since 



no Gg lactic foregrounds are known with 8 < —3 (sec Rc- 



ich & Reich (1988), who find that the low- frequ ency ra- 
dio s pectral index ne ar the NCP is ~ —2.7, and Banday 
& Wo lfcndalc (1991) , who present evidence that the syn- 
chrotron spectral index steepens to 8 ~ —3 at high fre- 
quencies). The integrated likelihood £(AT cm b) is shown 
in Figure [L7|. The rms obtained by maximizing £(AT cm b) 
is 



STrr 



82: 



9.1 



(37) 



where quoted errors define the 68% highest probability 
density (HPD) confidence interval. 




i:V' 



(uK) 



Fig. 17. — Integrated likelihood function £(AT cm |,) for the 
CMBR component. 



13. WINDOW FUNCTIONS 

The theoretical sky temperatures can be expanded in 
spherical harmonics T(xj) = J^i m a r t ^ m ( x i)' so that as- 
suming rotational symmetry, the expected value of the 2-pt 
correlation for fields i and j separated by an angle Xij on 
the sky is given by 



(38) 



where Pi are the associated Legendre polynomials, and 
Ci = (|aj"| 2 ) are the elements of the theoretical angular 
power spectrum (Peebles 1992). What we measure are the 
theoretical sky temperatures convolved with the telescope 
beam and switching strategy, i.e., AT^ = Ti * Bi (for the 
RING5M experiment, Bi(6,4>) is given by the function in 
Figure [l]) , so that defining Cy B = Bi eg) Bj, where Cg) de- 
notes cross-correlation, and letting * denote convolution, 
Eq. [58| gives 



I 



l)QPt(coBX«)*C, 



V B 



— Y(2l 



1)CiW{ 



(39) 



for the theoretical elements of the covariance matrix 
in Eq. ||. In CMBR parlance, the function = 
Pi (cos Xij ) * Cij B is known as a window function. For 
the RING5M experiment, a good approximation to the 
zero-lag (xu = 0) window function is give n by W" - 

exy\-l( l + 1)(t 2 ] (jSilk fc Wil- 



W*>m 



W aw 



where 



son 19801 , |Bond fc Efstathiou 1984|) , with a = 3'. 13, and 

| -2P,(cosx a ) + ^(cos2x fl ), with Xs = 22'.16. 



9= o 




Fig. 18.— Comparison of 0' (VK ; l,i ), 22' (W^ i+1 ) and 44' lag 

(W ; l,!+2 ) window functions for the RING5M experiment, character- 
izing correlations between neighboring fields. 



14. DATA CORRELATIONS 

14.1. Point Sources 

In general, point sources are not measured coinciden- 
tally with the Ring data, but are instead assumed to 
have constant flux densities and associated errors over 
timescales of a month (the typical time between VLA 
flux monitoring sessions during 1996) or longer. Since 
any source affects at least three fields through the double 
switching, subtraction of source contributions introduces 
correlated noise between neighboring fields, contributing 



(40) 



to the covariance matrix, where the bik are the beam 
weighting factors for the k th source in the i th field, and 
Ofc is the error associated with the flux density of the k th 
source. Although these covariances are included in the 
likelihood analysis, the effect is negligible, even for the 
1994-1995 RING5M data, for which source errors are en- 
larged to include variability. 

14.2. Noise Correlations 

An analysis of the long-term noise characteristics of the 
Ring5m data indicates the presence at 31.7 GHz of corre- 
lated noise between fields separated by 22'. The amplitude 
of this correlated noise is approximately 40 /iK, or 3—4% of 
the uncorrected noise in a single scan on a Ring field (see 
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§7.3). Analysis of subsets of the 31.7 GHz data suggest 
no obvious source for this component; we note however 
that its amplitude is of the same order as the observed 
season-to-season fluctuations in the 31.7 GHz mean levels 
suggesting that the same component may be 



7.4) 



(see 

responsible for both. 

While the origin of this correlated noise is not well un- 
derstood, independent confirmation of its presence can be 
seen through its effect on the data correlations; the mean 
Ring5m data set shows only half of the anti-correlation 
for nearest-neighbor fields expected for a 3-beam chopped 
experiment, Ci^±\jCn = —2/3. 

Our approach is to include the noise correlation as a free 
parameter in our model for the covariance matrix, i.e., we 
let 



= S 



V u ii 



(41) 



in Eq. |3q, where o~ n is the amplitude of the correlated noise, 
and a~~&re the variances in the field means, as before. As 
with the foreground spectral index, we integrate out the 
dependence on the noise correlation amplitude to obtain 
an estimate of the CMBR variance. Thus,the full form of 
the likelihood function £(AT cm b) in Eq. IS® is 



£(AT< 



cmb 



) = J J £(AT cmb 



P,a n )dl3da n . (42) 
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Fig. 19.— Integrated likelihood £(AT cmb ,<r n ) for the CMBR 
component extracted at the maximum likelihood value_qf /3. Con- 
tours are \ 2 = t T C~ 1 t for the 36 Ring5m fields (see §P-2|). Dashed 
contours are the 68% and 95% AT C111 |, — cr n confidence regions. Stip- 
pled region is the locus of 1-D 68% confidence intervals in AT CIn ^ 
at each value of cr n . 



In Figure |19|, we plot £(AT cmb , cr n ) for the CMBR com- 
ponent extracted at the maximum likelihood value of the 
foreground spectral index (j3 ~ —2.7), along with contours 
of constant x 2 , where x 2 = t T C _1 t (c.f. Eq. |34|). The 
peak of the likelihood function occurs for cr n = 41Tj l j l 'g /jK. 
Also shown is the locus of 68% confidence intervals on 
AT cm b which would be obtained if the correlated noise 
component were held fixed at the corresponding value on 
the y-axis. As can be seen from the figure, maximum like- 
lihood models which neglect these correlations are grossly 
discrepant with the data (x 2 ^> 34). 



15. BAND POWER 

The mass fluctuation power spectrum is often taken to 
be scale-invariant at small wave number, V(k) oc k (i.e., 
the Harrison-Zel'dovich spectrum), leading to a CMBR 
angular powe r spectrum w ith C[~ oc 1(1 + 1) on large an- 
gular scales ( Peebles 1992| ). As a result, the band power 

dBond 1995|) 

STi = y/l{l + l)Ci/2n (43) 

is expected to be flat at small and is most often the 
quantity predicted by theory. 

With the experimental rms obtained from the diago- 
nal elements of the covariance matrix in Eq. [39], ST Tms = 
VCwcmbi h can be seen that the band power is related to 
the rms by 



so that defining 



i ■ 



(44) 



(45) 



the weighted mean of the band power over the window 
function is given by 



5T, = 5T r 



(46) 



The integral of the window function is I(W") — 1.96, 
yielding 

M K, l e = 589±™l, (47) 



where l e — I(IW") / I(W") and errors on l e are the points 
at which the window function falls by e -0 5 . Errors on 
the band power are the 68% HPD confidence interval, and 
reflect sample variance, measurement error and 4.3% cali- 
bration error. 

16. DISCUSSION 

The goal of the RING5M experiment was to determine 
the CMBR anisotropy on 7' — 22' scales, corresponding to 
~ 10 Mpc at decoupling. We detect structure indepen- 
dently at 14.5 and 31.7 GHz well above the noise limits 
of the data. Observation of northern fields near transit 
ensures that the local environment is nearly identical from 
one field to the next; although the 14.5 GHz data show a 
parallactic angle dependence indicative of contamination 
by RFI, observations of the fields at both upper and lower 
culmination confirm that RFI contributes the same off- 
set to each field and thus does not affect our estimate of 
the variance. Numerous internal consistency checks (§0-|)) 
demonstrate that the structure observed in the RING5M is 
reproducible from year to year at both frequencies, while 
the strong correlation between frequencies confirms that 
the structure originates on the sky (§11). Careful charac- 
terization of both telescopes (§^|) and an extensive program 
of observations of calibrator sources (§|[) reduce our total 
calibration error to 4.3% — well below the sample variance 
of the experiment. 

With discrimination of foregrounds provided by our low 
frequency channel and accurate removal of point source 
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contamination using the VLA, we can state with confi- 
dence that 76% of the raw variance, or 88% of the source- 
subtracted variance at 31.7 GHz is due to the CMBR (see 
Table [?]). (A breakdown of the 31.7 GHz data set into 
CMBR, point source and foreground contributions is given 
in Table §) 

Our 14.5 GHz observations of the NCP have also re- 
sulted in the detection of an anomalous component of 
Galactic emission. The detection of this component is a 
cautionary tale for future CMBR experiments; whether 
the emission is due to high-temperature free-free emission 
correlated with Galactic dust (Leitch et al. 1997 ), or to 



some novel type of emission from the dust i tself (Ferrara 



& De;tmar 1994, Draine & Lazarian 19981), our results 



suggest that emission associated with Galactic IR cirrus is 
potentially a serious contaminant of small-scale anisotropy 
experiments even below 100 GHz. 

17. CONCLUSION 

A likelihood analysis of the RING5M data yields, for the 
CMBR component alone, a temperature rms of ST lms = 



82; 



-12.1 



UK, 



or equivalently, a band power of 5Ti = 
7o HPD confidence in- 



1+8.6 
-6.5 



+ l)Ci/2n = 59^ 

terval) at l e = 589^228- 

On 2' scales, the OVRO NCP 95% confidence upper 
limit of AT/T < 1.7 x 10~ 5 has recently been con- 
firmed by the SuZIE experiment, which set an upper 
limit of AT/T < 2.1 x 10~ 5 at nearly the same angu- 
lar scale ( |Church et al. 1997 ). While these results can- 
not strongly differentiate between varieties of closed uni- 
verses, as anisotropies are exponent ially damped at arc - 
minute scales for all of these models ( Bilk fc Wilson 1980| ) , 



the RING5M, NCP and SuZIE results together constitute 
a significant constraint on open universe models; models 
with h — 0.3 — 0.7 and fib = f7o = 0.2 over-predict small- 
scale power by 35% — 50% at ^ringsm, and by 10% — 35% 
over the NCP upper limit at Zncp- Open models with 
baryon density close to the lower bound allowed b y big 
bang nucl e osynthesis c alculations, Vl^h 2 ~ 0.015 (Wag- 



oner 1973, Dicus 1982), can reproduce the band powers 
observed at arcminute scales but severely under-predict 
the power observed by degree-scale experiments. 

Taken together, the RING5M and NCP results indicate 
a decrease in the angular power spectrum between I ~ 600 
and I ~ 2000. If the collection of data near I of a few 
hundred can be taken as evidence for a rise in the power 
spectrum toward small scales, then the RING5M band 
power is consistent with a peak in the power spectrum 
A x 2 fit to the data in Figure 19 for a range 



near I ~ 200. 

of model power spectra shows that the data are consistent 
with Qq = 1 in a A-model with H,t,h 2 = 0.015 and f^A = 0.7 
or a standard CDM scenario with Vl^h 2 = 0.0045, shown in 
Figure [lj]. Although increasing fit dramatically affects the 
amplitude of the first acoustic peak, the competing effects 
of damping at small scales and enhancement of compres- 
sional peaks with increasing means that the RING5M 
result cannot strongly constrain fib in a flat universe. 

Recent observations with the Cambridge Anisotropy 
Telescope (CAT) have resulted in a detection of anisotropy 
on angular scales directly comparable to tho se probed by 
the RING5M experiment flScott et al. 1996| ). The broad- 



with the CAT detection of ST h /T = 1.8+^ x KT 5 at 
I = 590. Given the result reported here, and the NCP and 
SuZIE results at higher I, there can be little doubt that 
the CMBR spectrum drops significantly between I ~ 600 
and I ~ 2000, as expected in flat cosmological models. 
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pointing out an error in the manuscript. We are grateful 
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the National Radio Astronomy Observatory, a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 



band power reported in this paper is in good agreement 



20 



CMBR Anisotropy 



Table 8 



Summary of Contributions to the 31.7 GHz Ring Field Means' 1 


Field 


ATraw (/iK) 


£raw (M K ) b 


ATsrc (mK) 


e src (/j,K) 


AT rm t (uK) 

' — CIIlD \f Ji ) 


frmh fuKI AT 


r nrr , (//Kl 
[Ore \r"^ / 


gf nr „ |//K) 

wore ) 


OV5M0024 


75.88 


17.63 


143.66 


17.07 


72.73 


50.34 


-140.52 


20.38 


OV5M0104 


-118.54 


17.73 


-146.49 


17.28 


-52.91 


45.35 


80.85 


19.14 


OV5M0144 


-43.42 


17.11 


-9.83 


0.84 


-20.42 


21.49 


-13.17 


5.43 


OV5M0224 


110.99 


15.77 


20.27 


1.15 


36.65 


20.09 


54.08 


5.36 


OV5M0304 


-86.21 


16.35 


-10.60 


0.90 


-55.68 


20.43 


-19.93 


4.86 


OV5M0344 


-19.82 


17.35 


1.72 


0.06 


-21.33 


21.64 


-0.22 


5.19 


OV5M0424 


90.33 


16.07 


10.75 


0.28 


114.66 


20.45 


-35.07 


5.71 


OV5M0504 


-13.56 


16.32 


1.92 


0.05 


-6.60 


20.38 


-8.89 


4.97 


(J V 0MU044 


— 11.28 


16.50 


2.62 


0.55 


—27.46 


20.55 


13.56 


4.82 


OV5M0624 


96.01 


15.60 


0.05 


1.87 


100.92 


19.72 


-4.96 


4.96 


OV5M0704 


132.13 


15.87 


5.30 


0.71 


174.83 


19.84 


-48.00 


4.77 


OV5M0744 


12.64 


15.29 


-1.23 


0.65 


18.66 


19.16 


-4.79 


4.75 


OV5M0824 


-122.33 


13.87 


1.80 


0.29 


-73.66 


17.33 


-50.47 


4.21 


OV5M0904 


145.15 


14.56 


4.62 


0.36 


99.82 


18.16 


40.72 


4.34 


OV5M0944 


45.37 


15.75 


-0.87 


0.35 


94.81 


19.68 


-48.58 


4.67 


OV5M1024 


33.42 


14.26 


13.68 


0.71 


5.21 


19.73 


14.53 


6.88 


OV5M1104 


0.56 


15.12 


8.38 


0.69 


-67.08 


19.19 


59.27 


4.93 


OV5M1144 


-78.24 


15.07 


-21.47 


1.17 


14.13 


20.43 


-70.89 


6.25 


OV5M1224 


109.97 


14.63 


58.15 


2.32 


-11.15 


23.51 


62.97 


9.28 


OV5M1304 


-62.25 


14.80 


-21.74 


2.38 


-38.47 


19.74 


-2.03 


5.81 


OV5M1344 


-83.76 


14.86 


18.59 


3.58 


-79.56 


19.29 


-22.80 


4.84 


OV5M1424 


-37.67 


15.06 


-28.83 


3.61 


-34.28 


19.93 


25.45 


5.97 


OV5M1504 


106.89 


15.74 


48.17 


6.14 


22.48 


22.14 


36.24 


7.48 


OV5M1544 


-46.12 


16.04 


-11.72 


2.61 


-27.78 


20.35 


-6.63 


5.07 


OV5M1624 


-68.42 


15.85 


-0.91 


0.29 


-71.48 


19.84 


3.96 


4.93 


OV5M1704 


-48.09 


15.69 


-9.81 


1.45 


-33.43 


19.70 


-4.84 


4.85 


OV5M1744 


35.44 


16.18 


22.33 


1.88 


-28.20 


20.50 


41.31 


5.26 


OV5M1824 


-126.52 


15.40 


-11.98 


1.70 


-102.90 


19.47 


-11.64 


5.03 


OV5M1904 


66.62 


16.01 


-1.79 


0.64 


60.59 


19.99 


7.82 


4.80 


OV5M1944 


-39.11 


15.28 


9.50 


0.53 


-32.86 


19.22 


-15.75 


4.94 


OV5M2024 


-192.57 


16.47 


-0.89 


0.71 


-169.52 


20.67 


-22.15 


5.17 


OV5M2104 


175.60 


16.14 


-0.97 


0.41 


157.75 


20.30 


18.83 


5.20 


OV5M2144 


-13.04 


16.62 


8.55 


1.72 


-59.85 


22.04 


38.26 


7.05 


OV5M2224 


176.05 


16.63 


12.56 


0.71 


149.43 


21.40 


14.07 


5.93 


OV5M2304 


-142.53 


17.61 


-7.96 


0.60 


-119.42 


21.99 


-15.16 


5.22 


OV5M2344 


-59.58 


18.60 


-105.50 


12.23 


11.36 


35.07 


34.57 


10.66 



a Equivalcnt R-J temperature. 

b Errors are la rms in the sample mean. 
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Fig. 20. — CMBR anisotropy measurements, shown with a A- model (solid line) with Q(,/i 2 = 0.015, Oa = 0.7, and Ho = 1, a standard 
CDM model (dashed line) with fl^h 2 = JL0045, and an open model (dot-dashed line) with fl),h 2 = 0.015 and Oo = 0.3. At bottom is the 
RING5M zero-lag window function (see §|13|), Indicated in bold are the OVRO RTNG5M detection of anisotropy at l e = 589 (this pappr) and 
the $5% confluence upper limit fmm tVi^nVRn Mrp oypprimnnt ( p^arllioarl gt a' ^QSQ > Other data paints are: COBE JHinaham et 

Tenerif P| (|HanCOck Ct a] 1W<), (pi.nHeT-s.Pn jt g .l™). Python jpiatt Ct al. 1922), ABm (Hp Rnrnardis et 




IRS ( |Ganga lgsg ), 

MAX (fl'anaka ct al. IQ Qffl ' Saskatoon 



were computed using CMHPS5T (scljak & Zaldarriaga 1996, icljak & ^aldarriaga 199Yr^aldarriaga et al. 199£ J 



et. al ma nsi i 



"CTiyl lgggS a5B Sn7IK ( 



Church ct al 1997| ) . Mode 
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